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SHARP ESTIMATES OF GREEN FUNCTION OF HYPERBOLIC BROWNIAN 

MOTION 

KAMIL BOGUS, TOMASZ BYCZKOWSKI, JACEK MALECKI 


Abstract. The main objective of the work is to provide sharp two-sided estimates of A-Green function 
of hyperbolic Brownian motion of a half-space. We strongly rely on recent results obtained by K. Bogus 
and J. Malecki [3], regarding precise estimates of the Bessel heat kernel of half-lines. 


1. Introduction 


The hyperbolic Brownian motion is the canonical diffusion on the hyperbolic space, with the half of 
the Laplace-Beltrami operator as its generator. It is related to many important objects, mainly to the 
Bessel processes and also to some functionals important in mathematical finance 0, m- Investigations 
of hyperbolic Brownian motion have long traditions; for a comprehensive survey, see e.g. m or 0,0, 
07 [Ml? [HI 7 [HI? [HI- Last years there was great interest in the process exiting a given domain. In the 
papers 0, m one can find a very useful representation of the Poisson kernel of a half-space, which led to 
provide its precise asymptotics. It turns out, however, that the most fundamental object of the potential 
theory is the A-Green function of the domain, A > 0. Unfortunately, its properties were not investigated 
until now in this context. The main objective of this work is to construct an adequate representation of 
the A-Green function of hyperbolic Brownian motion of the half-space {x = (ii,..., Xn) '■ Xn > a}, 
a > 0, and to provide its precise estimates. To achieve this, we employ Bessel processes > 0, 

of negative indices, killed at first hitting time of the point a. Recent paper [3] (see also [2]) contains very 
precise estimates for the transition density function pi ‘^\u, x, y) of the Bessel process BES^'^\ starting 
at the point x > a > 0 and killed at the first hitting time of a point a by the process. We briefly call the 
function pi '^\u,x,y) the Bessel heat kernel. 

Results of the paper are based on two important points: 

• the representation of in terms of the appropriate Bessel heat kernel; 

• sharp estimates of the above-mentioned Bessel heat kernel [3]. 


We first introduce some notation. If n denotes the dimension of the hyperbolic space, n > 
A > 0 we denote v = Furthermore, x = {x,Xn), y = {y,yn), Xn,yn > a, 

X = {xi,... Xn-i) G R””^. The main results of the paper are the following: 


Representation of G^ 


G^(x,y) = (XuPn) 


/i —u 


1 


/o (2-7™) 2 


— exp - 




2u 


pi ''\u,Xn,yn)du, 


where pi '^\u, Xn, y-n) is the appropriate Bessel heat kernel. 


2 and 
where 


Estimates of G^ 


For every p = 


n > 2, A > 0 and a > 0 we have 


G^{x,y) 


A,n / 2,Xnyn 


■-yV 


M-1/2 


1 A 


2{Xn - a){yn - a) 
\x - y|2 


1 A 


^Xjiyn 

\x - y\-‘ 


u-1/2 


whenever Xn^Vn > Here means that the ratio of the functions on the left and right-hand side is 
bounded from below and above by positive constants depending only on A and n. 
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Throughout the paper we rely on the representation of the joint density of an integral functional 
a\. = f* exp2(Bi ds and the geometric Brownian motion with drift exp(i?j = exp(_Bt — vt), 

given in m- An important point in our construction consists of application of the Lamperti Repre¬ 
sentation stating that the last coordinate of the hyperbolic Brownian, that is, the geometric Brownian 
motion exp(i?j~'^^), can be identified with the Bessel process with the time changed according 

to the functional 


2. Preliminaries 


In this section we collect some preliminary material. For more information on the modified Bessel func¬ 
tions we refer the Reader to [I] and nnj. For questions regarding Bessel processes, stochastic differential 
equations and one-dimensional diffusions we refer to El and to HH. 


2.1. Modified Bessel Functions. Various potential-theoretic objects appearing in the theory of Bessel 
processes, geometric Brownian motion or hyperbolic potential theory are expressed in terms of modified 
Bessel functions and K^. For convenience of the Reader we collect here basic information and 
properties of these functions used in the sequel. 

The modified Bessel function of the first kind is defined by (see, e.g. |1()| . 7.2.2 (12)): 

«-)-(2) £( 2 ) riuITlITT)' 


where d G R. The modified Bessel function of the third kind is defined by 


K^{z) = 


2 sin(d7r) 




R„(z) = hm R^(z) = (-1)"- 

ly—z 

We recall the asymptotic behavior of I^{z) at zero 


dl-ii 

dd 


dh 

dd 


i?=r; 


(see Uni, 7.2.2 (13) and (36)): 

( 2 ) 

, n G Z. (3) 


^ f(dTT)(l) ’ 

where g{r) ^ f(r) means that the ratio of g and / tends to 1, and at infinity 

Ii){z) = {2 'kzY^'^ e^W + 0{z~^)], z —>■ 00 . 


( 4 ) 

( 5 ) 


2.2. Bessel processes. The basic material concerning Bessel processes is taken from [T7], Ch. XL We 
begin with a definition of squared Bessel process BESQ^{x) started at cc > 0. It is defined as the unique 
strong solution of the equation 

dZ{t)=2^\Z{t)\d(3{t) + 5dt, Z{0) = x, (6) 

where I3{t) denotes a one-dimensional Brownian motion. Here d G R is called the dimension of BESQ^. 
For (5 > 0 the process Z{t) is non-negative and the square root in the equation (jS]) can be omitted. It is 
known that for 0 < d < 2 the point 0 is reflecting] for d = 0 it is absorbing. When d < 0 the situation gets 
more trickier: the process Z(t)^ when starting from a: > 0 attains the point 0 in finite time and, after 
that, becomes negative (and behaves like the process {—Z{t)) with positive dimension {—5) (see [H]). In 
this paper, however, we impose the killing condition at the point 0 and again call the resulting process 
Bessel process (of negative dimension). Thus, in our setting, the process Z{t) is always non-negative so 
we are able to take the square root. 

The square root of BESQ^(x^), a: > 0, is called the Bessel process of dimension <5 started at x and 
is denoted by {BES^{x)). We introduce also the index g = [5/2) — 1 of the corresponding process, and 
write BES^^^ (instead o BES^) if we want to use g (instead of 6). 

For g> 0 the probability density function of the BES^^'> (a:) semigroup is of the form 

p^^\t,x,y) = I (xy/t) for x > 0 (7) 

and 

p(M) )t, 0, y) = 2{2t)-f^-H-^^^+^'>T[g + 1)-^ (g) 

Note that the above-given formula describes the density considered with respect to Lebesgue measure. 
However, sometimes the symmetric version of the density can be more convenient. Then we have to 
deal with the reference measure m^^\dy) = m^^\y)dy, where We will switch to the 

symmetric version in Section [5.31 
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We denote by and the probability law and the corresponding expected value of a Bessel 
process BES^^\x) on the canonical paths space with the starting point a; > 0. Let Et = cr{i?s,s < t} 
be the hltration of the coordinate process Rt- Let us denote the first hitting time of the level a > 0 by 

Ta = inf{t > 0; = a}. 


We have the absolute continuity property for the laws of the Bessel processes with different indices 


dP 


(m) 




= — exp - 




ds 


, p 




- a.s. on {To > t}. 


(9) 


If 1 / > 0 then the condition (Tq > t} can be omitted. In this work we are interested primarily in the case 
of negative drift, which we denote by writing the index as (—/i), thus assuming in the sequel that /r > 0. 
Observe that the formula ([2]) implies that in the case of (—/i) the density function of BES^' X > Q 

is still of the form © with V = —/r. 


2.3. Exponential functionals. We now begin with a brief description of Matsumoto-Yor approach 
m, which is a starting point of our construction. The important point consists of the computation of 
the functional 


Hy) = E® 


exp - 


A* 


ds 

m 


\Rt = y 


To do this we apply the formula ([5]) to obtain 

X 


e1°) 


o-OiRt 


E(o) 


-oty 


© 


y. 

t \x 




exp - 


exp - 


x"^ + 




ds 

m 


2t 


\Rt 

h (y) ^(y)^y 




-ay y 


(!) 




t \X 


exp - 


x"^ + 
Yt 


'ImI 


(f) 


The above equality provides the desired formula 


^(y) = Ei°) 


exp I 


ds 

JEy 


\Rt = y 


Im (?) 

/o (y) 


Observe that the above formula defines the so-called Hartman-Watson distribution rjr with the parameter 
r = xy/t described in the terms of Laplace transform: 


fOO 

/ e~^ Vr{dt) = 

Jo 


Io{r) ' 


Below, we need a related density function, denoted by 9{r,t), defined as 

pOO 

/ = 7|^|(r) r > 0. 

^0 


( 10 ) 


3. Joint distribution of ^\b[ 


Let {Bt)t>o be the Brownian motion in R starting from x. We denote by P^ and E“ the corresponding 
probability law and expected value. Note that the starting point is now indicated in the superscript to 
distinguish it from the probability law and the expected value for the Bessel process with index y G R 
introduced previously. For y > 0 we denote by b[ = Bt — yt the Brownian motion with negative 
drift. We define 

pt 

4“'"^ = / exp(24”^^)ds. t > 0- 
Jo 

For y = 0 we will use the shortened notation At := Af'\ 

To find the joint distribution of (4 we proceed as follows. The crucial point is to determine 

the formula for the Green function G{x, y, o?Y = (di\ -f a^/2)“^, which is defined as a formal inverse, 
with El\ being the following Schrodiiiger operator 




iil 

2 dx^ 


{av. 


We first assume that y = 0. For x < y we obtain 


( 11 ) 


G(x,y,aV2) = 2/„(Ae©iL„(Ae©. 


(12) 
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It is easily seen that the Feynmann-Kac semigroup generated by the operator m is of the form 


pOO 

Jo 

ila for the transition density ' 
gx{t;x,y) = ^^\Bt = y - x] 


dy. 


lO \/2TTt 

This provides the following formula for the transition density of the above semigroup 

g-iv-xfl2t 


\/2Trt 


pOO 

/ e-^"/^^"^'^P°{AtGdu\Bt=y-x) 

Jo 


^-{y-xy/2t 

x/2'Kt 


Integrating with respect to t > 0 we find that 

pOO 

/ gxit;x,y)dt = G{x,y,a^/2) = 2 Ia{Xe^) Ka(Xe'^). 

Jo 

Taking into account the product formula for Bessel functions and the definition of the function 0(r, t) we 
obtain for x = Q 

'x^ey 


G(0,y,a72) = 


g-u/2g-AV2«g-A^ e^y/2'u 


g-AV2« g-1/2? g-e^G2« 


H/2 0 


-o^V2 0 - 


1 , du 
dt — 
u 

d^ 


= -« G2 


-AV2«g-(l+e^<')/2„W^ , 


C J ^ 


du 


pOO pOO 

Jo Jo 




'j2'Kt 


■dt. 


By comparing, we obtain for x = Q 


V\At€du\Bt = y) 


y/2'Kt 


= exp 


1 + e'^y \ „ f e.y \ du 


and for arbitrary x 


P%At€du,Bt = y) 


^-{y-x) I2t 

\/2-Kt 


= exp - 


^2x g2y 


2u 


u J u 


e^+y \ du 

—,0 —; 


(13) 


(14) 


which, after taking into account the drift formula (Girsanov) gives the following formula for the joint 
density of the variables {A[~^\b[~^'^) 


pa:(^(-G g du,B[->^'> G dy) = exp 9 ^ 


(15) 


4. Green eunction of { A ^~ y -\ exp { B ^~ y -'>)) 

We begin this section with providing the formula for the A-potential of the two-dimensional process 
(A(-M) , exp(i?( ^^)) together with its Laplace transform. The result can be found in the literature, but 
for completeness of the exposure and for convenience of the Reader we included it together with the 
proof in the following proposition. 

Proposition 1. X-potential of the process {A[ ^\exp{Bl is of the form 



where v = \j2X-\- pJ. The Laplace transform of is of the form 

>CQ^(a;,?/;-)(?'^/2) = - h{rx)K^{ry) = 2{xy)^ I^{rx)K^{ry)m'^~y‘\y). (17) 

y \yj 

m^' ^\y) = y ^ denotes here the density of the reference measure dm^ _ 
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^ e du, e dy)dt = exp (^) 


Proof. Taking into account the formula (TTSl) and (flUl) . the A-potential Q^{x,y;u) of {A\ ^\b[ is 
given by 

e=^ey\ 1 


u J u 


(18) 


where v = ^J2\ + y?. 

Using the formula for the product of modified Bessel function [see [S], p. 64 (37)] we can compute the 
Laplace transform of Q^(x, y; u) as a function of u 


CQ: 


poc poo 

\{x,y--){r^/2) = / & du, b[-^'> & dy)dt du 

Jo Jo 


M poo 


gxp 


2u 


e^ey\ 1 , 
— du 


u ) u 


a; \ M poo 


,-«/2 


exp - 


(re“)^ + (re^)^ \ fre^rey\ d^ 


2^ 




? J ? 


= 2 


\ M 


ey 


B{re^)K,{rey). 


Changing variables, we compute the A-potential Q^{x, y; u) of the process (A^ , exp(il( ^^)). We have 

for exp(_B(~^^) starting from x: 

px(^(-A.) g du,exp{Bl->^^) e dy) = ^ du,B[-^'^ € dlogy) 

V 




y \y 


exp 


+ y^^ „ fxy \ 1 


2u 


t] 

\ u / u 


and 


pOO 

Qt.ix,y,u) = / e~^*pr^{A[~^^ G du,exp{B[~^^) G dy)dt 

Jo 


I ( X 


= expi- 


-I- y ^ 
2u 


\ u J u 


The function Q^(-, •; w) is symmetric with respect to the reference measure dmf with the density 


d-M) 


(y) = 


1 


y 


2//+1 


When A = 0, we write Q^{x,y\u) instead of Q^{x,y\u). Thus, taking into account the previous formula 


we have 


Q%x,y\u) = 


pL — y 


Qv{x,y,u). 


The Laplace transform of u —> Q^ix, y, u) is given by 


CQ 


poo 

\{x,y,-){r^/2) = / e-^^^/^Q^{x,y,u)du 

Jo 




Jo y \y 

2 /x'' ^ 


exp - 


x^ + y^\ ^ (xy\ du 


2u 


r^y_a 
\ u J u 


B{rx)K„{ry) 

y \yj 

= 2{xy)^ I^{rx)K„{ry)m^~y-\y). 


□ 


To obtain formulas for the Green function we have to deal with the hitting time Tq. It is the first 
hitting time of a level a of the process exp(_B( ^^) starting from x > a. More precisely, we rather need 
to consider the density function qp°'{s) (where x is the starting point of exp(i?}~^^)) of the following 
stopped integral functional 

= f “exp(2B(-'^))ds. 
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We have 

where denotes the density investigated in the paper [7]. We have the following scaling properties 

of the density 

Considering an appropriate Schrodinger equation, we can write the Laplace transform of as follows: 




i=p(-y 4:'->)I = / exp(-L,),j(,)*. 


. 0 / K^{ra) 

Moreover, for every A > 0 we get in a similar way 




'[exp(-ATa)exp(-yx4^/>)] = (0 


x\^^ Ky{rx) 
a) K^{ra) ’ 


(19) 


where v = •^2A + p?. Further on we consider that the point a in the definition of q^'°‘(t) is fixed so we 
omit this superscript. Taking into account m we obtain for any positive Borel function /: 


E^ 


^[exp(-Ar,);/(47))] = E-[/(47))] = H f [s) ql{s) ds . 

\ Qj / \ Qj / J Q 


( 20 ) 


Theorem 1. The X-Green function of the of the process {A[ ^\exp(il( ^^)) killed at the first hitting 
time of the level a has the following form (for a < x < y) 


Gi,(x,y;u) = ( - ) [Q^{x,y]u) - Q^* ql{a,y]u)] 


with the corresponding Laplace transform 

f e"’'^“/^5;)(a:,?/;u)dM = 2 (xy)^ {[I^(rx)K^(ra) - K^{rx)I^{ra)) 

Jo K^(ra) 

= S^{rx,ra)m''~>"'>{y), 

where S^{a, (3) = I^ia) Kv{l3) — K„{a) 0 < /3 < a. 

Proof. We begin with the formula for the A-harmonic compensator of the process {A[ ^\exp{B^ ^^)), 
with respect to 


E"^[exp(-ATa); Q^(exp(B(_^ ^)), j/; + u)] = 


E"'[exp(-Ara); Q^{a,y; + u)] 


/ a 


© 


poo 


Q„{a,y;s + u) ql{s)ds 


Qy * q(iia,y;u). 


Next we compute the Laplace transform of the compensator with respect to the variable u 
CQ^*ql{a,y]-){r'^/2) = TQ^{a,y\-){r^(2) ■ Tqlfr'^(2) 


= -{-] I, 

y \y 


2 fx 

y \y 


{ra)K^{ry) ■ (^0 
K^{rx)K^{ry) 


= 2{xy)'' K„{rx)K„{ry) 


x\'^ Ki,(rx) 
aJ Ky{ra) 
Iy{ra) 

K„{ra) 

Ki,{ra) 


©)• 
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The A-Green function of the process {A[ ^\exjp{Bl killed at the first hitting time of the level a has 
the following form (for a < x < y) 


/ x \ 

G^{x,y;u) = f-j [Q^{x,y;u) - Q^*q^{a,y;u)] 


poo 


yj Jo 


[Q,,{x,y;u) - Qi,{a,y;u- s)l[o,„)(s)]g;^(s)ds. 


Its Laplace transform is 


[ e '^^^/'^g^{x,y,u)du = - f-) I^{rx)K^{ry) - - (K^{rx)K^{ry) (21) 

Jo y \yj 


y \y 


K^{ra) 


2 {xyY {{Iu{rx)K^{ra) - K^{rx)I^{ra)) ^\y) (22) 

K„(ra) 

2 {xyT Sv{rx, ra) (y), 


with S^{a,P) = I^{a) K^iP) - K^{a) I^{I3), 0 < P < a. 

The above-given result immediately implies the following 

Corollary 1. 


(23) 

□ 


G^ix,y;u) 


Gv{x,y]u). 


5. Green function estimates 


5.1. Hyperbolic Brownian motion in H". We consider the half-space model of the n-dimensional 
real hyperbolic space 

H" = {a; = (xi,..., x„) G R" : Xn > 0}, n>2, 


with the Riemannian metric 


, 2 dxl -I- ... -I- dXn_i + dx„ 

ds = -^- 


The metric induces the hyperbolic distance on H" described by 

I |2 

cosh(dH"(a;, j/)) = 1-f 

and the corresponding canonical (hyperbolic) volume element 

dxi ... dxn-idxn 


x,y € 


(24) 


dVn = 

Xn 

where dxi ... dxn-idxn denotes the Lebesgue measure in R”. The Laplace-Beltrami operator takes then 
the following form 

” ^2 Q 

Ah- = x" ^ ^ - (n - 2)a;„ —. 


fc=i 


We define the hyperbolic Brownian motion (HBM) X{t) = {Xi{t),..., Xn{t)) as a diffusion on H" 
with the generator ^Ah". Moreover, we introduce the hyperbolic Brownian motion with drift, i.e. the 
diffusion X''(t) = {X ^,..., X'^{t)) on H" having the half of the operator 

" a2 Q 

as its generator, where u > 0. Notice that 

11 d 

= -Ah« - (u - (n - l)/2)a;„^^, 

and consequently, for y, = (n — l)/2 we have A^ = Ah" and we go back to HBM, i.e. X = X^. We will 
fix the notation y = (n — l)/2 in this section. 
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5.2. Representations of the hyperbolic Green fnnction of the horocycle. Our main objective is 
to study the properties of the Green function and the A-Green function of a half-space (or equivalently 
the interior of the horocycle), i.e. the set 

D = {x & H" : Xn > a}, 

where a > 0. We denote by the first exit time of X'' from D 

Ti, = inf{t > 0 : X'"{t) ^ D} — inf{t > 0 : X'^{t) = a}. 

Since the last coordinate of the hyperbolic Brownian motion with drift has the same law as the cor¬ 
responding geometric Brownian motion, we can easily deduce that is finite almost surely whenever 
a > 0. For every A > 0 we define the A-Green function G^{x, y) of D for the hyperbolic Brownian motion 
with drift u as the integral kernel of the Green operator 

G^/(^) = r < Ta, fiX'^m dt= [ f{y)Gt{x, y)dVn{y), 

Jo J Da, 

for every x,y G D and any Borel function /, which is non-negative or bounded. We do not indicate the 
value of a in the notation of the corresponding objects, i.e. a is assumed to be fixed positive number. 
We also recall the reference measure formula dVn{y) = yn'^dydyn = m^~^\dyn)dy, where y = (n— l)/2. 
The A-potential operator and the A-potential kernel U^{x,y) are defined by 

Ul^f{x)= r e-^*E[f{X’^{t))]dt= [ f{y)U^{x,y)dVn{y), x,yGW^. 

Jo J Da 

For A = 0 we obtain the Green function of D and the potential kernel, which will be simply denoted 
by G^{x,y) and U^{x,y) respectively. Finally, ii v = y = (n — l)/2 we will omit the subscript v in the 
notation, i.e. G^{x, y) and G{x, y) are the A-Green function and the Green function of the set D for the 
hyperbolic Brownian motion. 

The hyperbolic Brownian motion with drift can be represented by the classical Brownian motion 
with time changed by an integral functional of geometric Brownian motion and the geometric Brownian 
motion itself. More precisely, there exists the standard Brownian motion B = = {B{t),Bn{t)) in 

R" starting from x = (i, loga;„), Xn > 0 such that 

X-'it) = {X''{tlXl{t)) = 


where si — vt, A[ = f*exp(2Bi ‘'\s))ds. Since the processes B(t) and the two- 

dimensional process {A\. *^\exp(i?i are independent and the first exit time Ti, depends only on 

the last coordinate of the process, the above-given representation implies the relation between hyperbolic 
Green function of the horocycle and the A-Green function of , exp(i?i~^^ (t))) introduced in Section 

SI More precisely, the A-potential operator for the hyperbolic Brownian motion X and non-negative 
function / can be written as 

pOO 

U^fix) = / 

nOC nOO 

/ / Wf{B{w),y,,)P^-{A[-^'^ Gdw,exp{Bl>^\t))Gdy^)dt. 

Jo Jo 

Since the integrand is non-negative, we can change the order of integration to get 
poo p poo exof_ 

U^f{x)= / / f{y,yn)— - , G dw,exp{Bi-^^\t)) G dyn)dtdy 

Jo Jo ( 2 ttw) 2 Jo 

= [ [ Q^ixn,yn;w)dwdy. 

JR^ Jo (ZTTW) 2 

Consequently, using the relation 




Qi^iXn,yn',w) 


and the formula for the hyperbolic volume element we obtain the corresponding formula for the A- 
potential kernel 


U^{x,y) = y 


n 


fj — u pOO 


{2'!Tw)^~ 


'OO 


Qv{Xn,yn]w)dw. 


(25) 
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Similar computation gives the formula for the A-harmonic compensator 

— U 




a 


Vr 


= 


= ,,2m+1 

Un 




= I - 

\y. 


fi — y 




^oo noo 


exp(- 




\B{s)-y\^ 


-At.®^P( 2w ) 


/O 


{2ttw)~ 


Qu{a,yn\w) dw 


0 ^0 

roo .oo . exp(-H^) exp(- 


/ 


, , n — 1 

{2'kw) 2 

2w 


— - -Qi^{a, yn] w)ql^ (s) dw ds 




0 Jo {2 ttw) 2 


-1 


(27rs)' 


^ 71 — Qj/(a, ?/„; (s) dwdsdz 


y-—y poo poo 


exp(- 




) 


, , , , Qy{a,yn]w)ql''{s)dwds 

Jo Jo {2tt{w + s)) 2 

/■°° exol- 


Jo (2?™) 2 

The above-given relations enable us to provide the following representation formulas for the A-Green 
function G^{x,y) of a half-space for the hyperbolic Brownian motion: 

Theorem 2. The X-Green function of D is given by 


G^{x,y) = yl^^^ — - ^-1 G^{xn,yn;w)dw, 

Jo (2™) — 

where Q^(xn,yn]w) is the Green function for the process {A\ ^\exp(_B( ^^))- 


(26) 


The above-given formula together with the relation between a Bessel process and a geometric Brownian 
motion established by Lamperti relation lead to the following result, which will be crucial in estimating 
the A-Green function. 


Theorem 3. 


G^{x,y) = {XnynY 


— exp 


'-y? 

2u 


pi ’'\u,Xn,yn)du, 


(27) 


Jo (2?™) 2 

wherepi (u, a;^, y„) is the density function (with respect to the speed measure'mJ~'^\dyn) = dyn) 

of the transition probability of the Bessel process BES^~G [xn) killed at the first hitting time Ta- 

Proof. From the formula and the fact that 


we obtain 


Qyn: — i ) Qui.Xn^yn^'xJ) 
. yn 


Gix,y)=[ — ) Gi,{x,y) x,y € D, 
. yn 


where u = ^y2X + and y = {n — l)/2. To prove our theorem, we use the following fact: after changing 
variables in the geometric Brownian motion according to the formula t = Uu where Uu = inf{s > 
0; Ai~G > u} we obtain that exp(Bi~G'j = BES^~G is a Bessel process with index (—u), with Ai G 
before. Applying this relation we obtain that Ta = inf{s > 0; exp(i?i~'^^) = a} is the first hitting time 
of the Bessel process BES^~G of the level a and ■ Applying these relations, we obtain for a 

positive Borel function /: 


GJix) = E" / l^t<r.yfiB^,-.yexp{Bi-G{t)))dt 
Jo * 

pOO 1 


^OO pOO 


E''[fiBu,yn)]pi G{u,Xn,yn) 


dyn 


10 Jo 


,,2,,2:^-1 
duPn 


du 


/R" Jo {2ttu) 2 


^ — exp (- ] p[_ G (u, X„, yn) yl^ fiy,yn) dudVniy) 


2u 
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Collecting all together provides the result. 


□ 


5.3. Sharp estimates of Green functions. We begin with providing the two-sided uniform estimates 
of the A-potential kernels. 


Proposition 2. For n > 3 we have 


U\x,y) 


\x - yp 


m-1/2 


1 A 


‘^XyiVu 

|a; - yP 


v+l/2 


x,y G 


where fi = and v = y?. Moreover, if n = 2 we have 


U\x,y)^{ lA 


2^Xjiyn 

\x - yp 


v+l/2 


, X,y& 


Proof. Using the formula (PSI) together with the asymptotic description of Q^{xn,yn',t) given in ([I]) and 
([^ and the fact that y = we get 


U^(x,y) « (x„y„)^ du + (x„y„) 


\I1 + V 


1 

e 2“ du. 


Jo u 

Making the substitution = s in both integrals, we arrive at 

Using the following asymptotic description of incomplete gamma functions (see Lemma 12 in [5j) 


pb 

/ s“e-*ds Ri (1 A6)“+i, b> 

Jo 


s“e ®(is r: (a-I-l)“e a > 0, 


we obtain that for P > 1 that 

2x„y„ — 


U\x,y) 


2x„y„ 


\x - yp 


e 2a:nl/n -|- 


^^nUn 

|x-yp 




‘^^nUn 

|x-yp 


fl+U 


Finally, whenever < 1 we have 


U\x,y) 


‘2Xnyn 

|x-yp 




-f 1. 


Note that if n > 3 then the potential behaves as the first summand and in the case n = 2 as the other 


one. This ends the proof. 

Theorem 4. For every y = n > 2, A > 0 and a> 0 we have 


□ 


G\x,y) 


A.n / 2x„y„ A'" ^/^ / 2(x„ - a)(y„ - a) 


1 A 


1 A 


‘2‘Xnyn 
\x - yp 


v-\/2 


|x-yp/ V |a:-yp 

where x = (i,x„), y = (y,y„), x„,y„ > a and v = ^J2\ F y^. 

Proof. The scaling property of (HBM) gives 

G^{ax,ay) = G^{x,y), x„,y„>l 

where Gg„ denotes the A-Green function for the set Di. Consequently, we will further consider a = 1 
and omit the sign ~. We recall the recent result giving the sharp estimates of the transition density of 
the Bessel process killed when leaving the half-line (a, oo), i.e. it was shown in [3| that for every u 7 ^ 0 
we have 


(—1/)/. I U 1 A l)(y" Pi A ^riyn\'' 1/^ / , 

Pa (t:Xn,yn)'^ \ ^ t /V^t/ (Xnyn) 


v-1/2 


1 

Vt 


exp - 


(Xn - y-af 


2t 


, (28) 
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whenever t > 0 and XmlJn > 1- Note that in [3] the density was considered with respect to the Lebesgue 
measure and here the reference measure is m{dyn) = dyn- Splitting the integral (1371) representing 


G^{x,y) into three parts we get 


G\x,y) = 


rXnV 


pXnVn rO 

/ / 

J (Xr^—l)(Vr,—l) Jxr, 


^(rCn-l)(yn-l) '^XnVn 

~ Ji{x,y) + J 2 {x,y) + J3{x,y), 

where, by substituting w = \x — yp/(2t) we get 


{XnVnY '' 

(27rt)("-i)/2 


exp - 


''\xn,yn)dt 


Ji{x,y) = {XnynY [ 

Jo 


(rCn-l)(yn-l) 


‘J^nVn 


i/2-l 


I \x-y\^ 

2(x„-1)(b„-1) 


and similarly 


y) = (' 2 (a^n - l){yn - 1) /•^Cn-uo„-u 

\\x-y?) \x-y? 


t-"/^exp(-^)dt 
wr^/2-2e-^dw 

-J/l" 


(29) 


(30) 


U-i/U 


Note that if \x — y\'^ < 2(x„ — l){yn — 1) then the integral in (1331) behaves like a constant and 


Ji{x,y) 


‘^Xfiyn 


n/2-1 


lx - y|2 

Moreover, since |x — j/p < 2(x„ — l){yn — 1) < 2xnyn we have 

2(x„-1)(2/„-1)-^ 


J2{x,y) 


2xnyn\"'^ ^ 2(x„-!)(?/„-1) /■2(x„-l)(„„-l) 


|x-y|^ 

2Xnyn 

|x - J/p 


and since 


n/2-l 


( 2XnynY""'' ^ 2(x„ -!)(?/„- 1) /■2x„;„ 

— Jo " 


\x - yp 
|x-yp 


2(x„ - l)(y„ - 1) 


ra/2-l 


< 


2x V ^ 

^^nyn 


■-y|' 


< 


< 


\x - yp 
{Xn - l)(yn - 1) 


^riVn 


^^nUn 

\x - y|^ 


^/2-l 


\x-y\" 


2{Xn - l)(yn - 1) 


i/2-l 


< 


2Xnyn 

|x-y|^ 


1 / 2-1 


which ends the proof in this case. 

In the second case, i.e when 2(x„ —l)(y„ —1) < |x —yp < x„y„, the integral part of (l30l) is comparable 
with a constant and consequently, J 2 {x,y) dominates the other parts. To see that note the asymptotic 
formula for the incomplete gamma function 


w°' '"dw ~ z“ z —>■ 00 . 


a—1 ^—z 


(32) 


which gives that the integral part of (USD decays exponentially as |x — yp/(2(x„ — l)(yn — 1)) grows to 
inhnity and consequently Ji{x,y)/J^ix^y) is bounded from above in the considered region. Moreover, 
the above given estimates of J 3 {x,y) and J 2 (x,y) give 

T ( \ ^ (^" ~ OiVn ~ 0 ^ ^ j f \ 

J3{x,y) < -< 1 « J 2 (x,y). 


^nUn 
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The final part of the proof relates to the case when 2{xn — l)(2/n — 1) < |a; — 2/P and XnVn < |a; — yp. 
Then, by the other assumption, the integral part of (1511) behaves like a constant and to finish the proof 
it is enough to show that J 3 {x, y) dominates Ji(x, y) and J 2 (x, y) in that case. Indeed, using (1!?^ we get 


J2[x,y) < 


‘^XnVn \ ^ 2(Xn - l)(y„ - 1) 

k - 2/P 


< 


\x - J/p 
‘2‘Xnyn 

k-2/P 

‘2Xnyn 

\x - 2/P 


' \x-y\^ 
"^^nyri 




2{Xn - l)(2/n - 1) (\x-y\ 


\x-y? 


‘2Xnyn 


exp - 


/i+jy-i 


2{Xn - l)(2/n - 1) f\x-y\ 


|a; - 2 /P V ‘^XnU. 
where ci(n) = w^~‘^e~^dw)~^ • sup,j,>i /2 Finally, 


exp - 


\x - 2 /P 

2'Xnyn 

jx-y? 

‘2‘XnyrL 


< ci(n)J3(a:,2/), 


Jiix,y) < 


< 


< 


2‘Xjiyn 

\x - 2 /p 

2^Xnyn 

\x - 2 /p 

2^Xnyn 


i/2-l 


\x - 2 /P 


2{Xn - l)(2/„ - 1) 


n/2-2 


exp - 


|a:- 2 /P 


n-2 


2{Xn - l){yn - 1 ) 


|a;-2/P 


\x - 2/P 

2(x„-1)(2/„-1) 


2{Xn - l)(2/n - 1) 


2(x„ - l)(2/n - 1) 
n—2 

exp 


-2/P 


I®-2/P 


-2/P 


2{Xn - l)iyn - 1) 


2{Xn - l)(2/„ - 1) 

n — 2 / I 12 

\x-yr 

2{Xn - l)(2/n - 1) 


exp 


\x - 2/p. 

< ci{n)J3{x,y). 

This ends the proof. □ 

The above-given result can be stated in terms of the hyperbolic distance in the following way. 
Corollary 2. For every n > 2 and X > 0 we have 


\ , . 'Xd'‘ 

G\x,y) ^ 


1 


sinh^^ ^ (c/h" (x, y)/2) cosh‘'+^/^((iHn {x, y)) 


1 A 


(1 A <Ia(x))(l A 5a{y)) 
1 A (x, y) 


, x,y G D, 


where Sa(x) denotes the (hyperbolic) distance of x to the boundary of set D, or equivalently 


A A<5a(2/))A ..A 


G^(x, 2/) 1 A „ 

V 1 A (x, y) 

Proof. First observe that by we simply have 

|2\ -1 


u (x,2/), x,y G D. 


1 A 


‘2‘Xjiyn 


\x - 2/P 


1-b 


X - 2/P 

‘2‘Xnyn 


= cosh \(iH"(x,2/)). 


Moreover, (1^^ also implies that 

‘2‘Xnyy 


-1 1 -2 

^ _ ^|2 = (cosh(dH-(x, 2/)) - 1) =-sinh {dn'^{x,y)/ 2 ). 


Since the lines perpendicular to the boundary of the set D are geodesis of the hyperbolic space, the 
distance of x = (x,x„) to the boundary of D is realized as a distance between x and (x, a). Thus we 
have 


cosh 5 a (x) = 1 -b 


{Xn O,) 1 f Xn L , X ^ \ 1 / / ^ 


Consequently 


a 


— 1 = — 1 ~ sinh(5a(x), — = ~ cosh 5a(x). 

a 


Using this we can write 

(x„ - a)(2/n - a) (x„/a - l){yn/a- 1) 2x„2/„ sinh5a(x) sinh5a(2/) 


■-2/P 


2 x„ 2 /n/c 


|x - 2/p cosh5a(a:) cosh5a(2/) 2 sinh^(dH" (a;, 2/)/2)' 


To finish the proof note that tanh5a(x) ~ 1 A 5a(x) and 
(1 A 5a(x))(l A Sa{y)) 


1 A 


sinh^(dH"(a;, 2 /)/ 2 ) 


COShdHn(a,^y) Ri (1 A 5a(x))(l A 5a(2/)) 



































SHARP ESTIMATES OF GREEN FUNCTION OF HYPERBOLIC BROWNIAN MOTION 


13 


whenever dm" ix,y) > 1 and 


(1 A i5a(a:))(l A i5a(y)) \ 
sinh^((iH»(a;,i/)/2) J 


cosh d^rt (^x,y) 


(1 A (^a(x))(l A Sqjy)) 
d^„ (x, y) 


if only dM'^{x,y) < 1. Collecting all together we get the result. 


□ 


We end with the following conjecture stating that the form of the above-given estimates should be 
the same for more general smooth subsets of H". Obviously, the conjecture is true if we assume that the 
considered set D is bounded in hyperbolic metric. It follows from the fact that the potential theory on 
bounded subsets of H" is comparable to Euclidean potential theory, since the Laplace-Beltrami operator 
is strongly elliptic operator on such domains. 


Conjecture 1. For every domain D C H" we have 


GD{x,y) 


n,D 


1 A 


(1 A (5a(a:))(l A da{y)) 
1 A (x, y) 


U{x,y), 


x,y G D. 
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